Hybridization and spin decoherence in heavy-hole quantum dots 
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We theoretically investigate the spin dynamics of a heavy hole confined to an unstrained III-V 
semiconductor quantum dot and interacting with a narrowed nuclear-spin bath. We show that band 
hybridization leads to an exponential decay of hole-spin superpositions due to hyperfine-mediated 
nuclear pair flips, and that the accordant single-hole-spin decoherence time Ti can be tuned over 
many orders of magnitude by changing external parameters. In particular, we show that, under 
experimentally accessible conditions, it is possible to suppress hyperfine-mediated nuclear-pair-flip 
processes so strongly that hole-spin quantum dots may be operated beyond the 'ultimate limitation' 
set by the hyperfine interaction which is present in other spin-qubit candidate systems. 
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Heavy holes (HH) confined to semiconductor quantum 
dots (QD) have attracted rapidly growing attention over 
the last years for their potential applicability in spintron- 
ics and as qubits for quantum information processing de- 
vices. The spin states of confined HHs feature very long 
spin relaxation times [IH3] and are suspected to be ro- 
bust against spin decoherence. Typically, for spin qubits 
operated at sub-Kelvin temperatures, the main source of 
decoherence is the interaction with the nuclear spins re- 
siding in the host material and the inhomogeneous broad- 
ening of the nuclear magnetic field (Overhauser field) @]. 
For HHs, the form of the nuclear-spin interaction is pre- 
dominantly Ising-like [5], in contrast to the Heisenberg- 
type interaction of electrons. Hole-spin QDs in p-type 
GaAs/AlGaAs heterostructures have already been real- 
ized experimentally and operated in the few-hole regime 
[BHH]. Experiments in self-assembled InGaAs quantum 
dots have shown the possibility to initialize and read out 
the spin state of a HH with high fidelity [9] , and ensemble- 
spin decoherence times T% on the order of hundreds of 
nanoseconds have been measured |10j . 

Several possibilities to suppress decoherence due to 
inhomogeneous broadening have been proposed, one of 
which is to prepare the nuclear spins in a so-called nar- 
rowed or frequency- focused state [TTl - fT5] , where the bath 
is prepared in an eigenstate of the Overhauser operator 
(see text below Eq. ([6])). On the experimental side, enor- 
mous progress has been achieved in preparing such nar- 
rowed states [HHS], which have been shown to persist 
over astonishingly long timescales exceeding hours |15j . 
For electrons interacting with a narrowed nuclear bath, 
spin decoherence happens due to nuclear pair flip pro- 
cesses induced by the transverse hyperfine interaction, 
and the associated single-spin decoherence time T 2 can 
be several orders of magnitude longer than the ensemble- 
spin decoherence time T| [251 H3]- F° r HHs, with their 
predominantly Ising-like coupling to nuclear spins 5J, 
this transverse interaction (perpendicular to the Ising 
axis) can be expected to be very small, potentially lead- 
ing to very long single- hole-spin decoherence times T 2 . 

In this Letter, we study the spin dynamics of a HH 



confined to a III-V semiconductor QD and interacting 
with a narrowed nuclear-spin bath. We show that band 
hybridization leads to non-Ising (transverse) terms in the 
hyperfine Hamiltonian, whose magnitude depends on the 
geometry of the QD. This transverse coupling induces 
nuclear pair-flip processes, leading to fluctuations of the 
Overhauser field and to exponential single-hole-spin de- 
coherence. We show that for typical unstrained quantum 
dots the associated timescale T 2 has a lower bound on the 
order of tens of microseconds and that it can be tuned 
over many orders of magnitude by changing external pa- 
rameters such as the applied magnetic field. Thus, it is in 
principle possible to operate hole-spin QDs in a regime 
where the hyperfine interaction is practically switched 
off and where other decoherence mechanisms, such as 
nuclear dipole or spin-orbit interactions, will become rel- 
evant and, hence, experimentally observable. 

We start from the 8x8 Kane Hamiltonian describing 
states in the conduction band (CB), heavy-hole (HH), 
light-hole (LH) and split-off (SO) bands of bulk III-V 
semiconductors (see Appendix C of Ref. [H]). The Kane 
Hamiltonian can be 'folded down' to an effective 2x2 
Hamiltonian whose eigenstates describe the spin states 
in the band of interest, where the admixture of neigh- 
boring bands is taken into account perturbatively |24j . 
Using this procedure, we find the following hybridized 
HH pseudospin states (see Appendix |A| : 

|*±)~ A/(|4 h ;0h°H>I±HH> 

T A CB |4 B ; 0±cb) ± Alh |wl H ;<0±lh))- (1) 

Here, we have assumed a parabolic confinement poten- 
tial defining a QD with lateral and perpendicular con- 
finement lengths L and a z , respectively, and Af enforces 
proper normalization of the wavefunctions. The condi- 
tion for the validity of Eq. ([!]) is given by a z -C L, 
which is needed for the perturbation expansion on the 
Kane Hamiltonian. The amount of CB and LH ad- 
mixture is determined by Acb = i/3cBP/V2LE g and 
Alh = \/3f3- L Yi-f3a z L/2y/2j2(L 2 -a 2 z ), respectively, where 
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P is the interband momentum, E g is the band gap, 72,3 
are Luttinger parameters, and /3cb, Pixl account for the 
difference in effective masses between the bands (see Ap- 
pendix [A| . 

Near the T-point, the spin-orbit-coupled states can be 
approximated by |u^ B )|± CB ) - \s) \ |mhh)I ± hh> ^ 

|p±)lt4), N 



LH, 
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lh; 



(V2\ Pz )\t,i)T\p±)\U))/V3, 
in terms of s- and p-symmetric Bloch states (\p±) = 
\p x ) ± i\p y )) and real-spin states | t>4-) with respect 
to the growth axis [23] • The envelope functions ap- 
pearing in Eq. ([I]) are defined via their position rep- 
resentations (r|$jf) = 4> % ^-{z)(0{x,y) (i = 0,1, j = 
0,±), where <$ l (x,y) = (j)° a (x)<j) a (y), (f>t ll (x,y) = 
{4>i{x)f a {y) ± i<t>° a (x)f a (y))/V2, and <fc(x) is the n th 
harmonic-oscillator eigenf unction in band a. Due to 
terms appearing in the Kane Hamiltonian which are lin- 
ear in the crystal momentum k and which couple neigh- 
boring bands, the admixture of CB and LH states fea- 
tures excited-state envelope functions. This has profound 
physical consequences which will be discussed below. The 
split-off-band contribution to the HH states is very small 
and has thus been neglected in Eq. ([lj. 

There are three interactions that couple an electron (or 
HH) to the spins of the surrounding nuclei: the Fermi 
contact interaction h^, the anisotropic hyperfine interac- 
tion h\, and the coupling of orbital angular momentum 
to the nuclear spins /I3, which read (setting H = 1) [2"5] : 
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hS = 
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Here, js = 2mbi lj k = 9j k ^Ni Ms is the Bohr magneton, 
gj k is the nuclear g-factor of isotopic species jk at lattice 
site fc, fj,N is the nuclear magneton, = r — is the 
electron-spin position operator relative to the fc th nucleus 
with spin If., d ~ Z x 1.5 x 10~ 15 m, Z is the charge of 
the nucleus, and n& = r^/r^. S and L^ = r*, X p denote 
the spin (m§ — ±1/2) and orbital angular-momentum 
operators of the electron, respectively. 

In order to derive an effective spin Hamiltonian for the 
HH, we take matrix elements (* T |^i + ^2 + ^I^V) = 
H TT ' (t, t' = ±) with respect to the hybridized HH wave- 
functions |T]). Due to the (5-function in Eq. |2|, only the 
CB admixture contributes to the Fermi contact interac- 
tion, since p-states vanish at the positions Rft of the nu- 
clei. On the other hand, the terms in Eq. ([!]) associated 
with HH and LH states contribute to matrix elements of 
Eqs. ([3]) and Q, while the CB admixture does not con- 
tribute due to symmetry (h^) and vanishing orbital an- 
gular momentum (ft.3). Adding up all contributions, and 
taking into account a Zeeman term due to a magnetic 
field B along the z-direction, we find the following ef- 
fective spin Hamiltonian describing the hole-nuclear-spin 



interactions: 



H = {b + h z )S z + -(h+S- +h-S + ). (5) 

Here, b — ghP^sB is the Zeeman energy of the HH, 
cjh — 2 is the HH g-factor along the magnetic-field di- 
rection z, fiB is the Bohr magneton, and S is the HH 
pscudospin-1/2 operator. The Overhauser-field compo- 
nents are defined by h z — J2k ^Jk an d ^ = J2k ^k^k 
(/jf = If ± Hi), wnere A-% and A k denote the longitu- 
dinal and transverse hyperfine coupling of the HH to the 
fc th nuclear spin, respectively. The flip-flop terms in Eq. 
([5]) couple the HH pseudospin (±3/2) states through ad- 
mixture with CB and LH pseudospin (±1/2) states, such 
that flip-flop processes with 1=1/2 nuclear spins pre- 
serve the total angular momentum. 

The hybridized states in Eq. (lj are predominantly 
HH-likc. In Ref. [5] it has been shown that taking ma- 
trix elements of the Hamiltonians with respect 
to pure HH states (i.e., neglecting band hybridization) 
results in an Ising Hamiltonian h z S z . The longitudi- 
nal coupling constants are thus dominated by the HH 



contribution, At 



A z k HH , and the transverse (lion- 



ising) terms in Eq. ^ are only due to hybridization 



with CB and LH states, At 
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plicitly, the longitudinal and transverse coupling con- 



stants are given by A 



fc,HH 
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and 



A ± 



k,GB 

- A ik v o\^i(zk)\ 2 <l>^(xk,yk)(l)±(xk,yk), respec- 
tively, where vq is the volume occupied by one nucleus 
and A£ is the hyperfine coupling strength of isotope 
jk associated with band a. Introducing the average 
A a = J^j v jA J a , where vj denotes the abundance of iso- 
tope j, we estimate Ahh — — 13peV [5], Aqb — 0.15/xeV, 
and Alh ~ 0.05/ieV for a GaAs QD with L = lOnm 
and a z = 2nm. In contrast to the interaction of an elec- 
tron with nuclear spins, the hole-nuclear-spin interaction 
given in Eq. ^ is highly anisotropic. 

We now study the dynamics of the transverse spin 
component S + describing the coherence of the HH pseu- 
dospin states. To this end, we use the Nakajima-Zwanzig 
master equation [1 lj 



(S+) t =iio n (S + ) t -t dt^{t-t'){S + ) r , (6) 



where uj = b + h z , w\n) = uj n \n), and \n) denotes a nar- 
rowed state of the nuclear-spin system (note that for a 
non-narrowed bath, the HH decoherence would be dom- 
inated by the Ising part of Eq. ([5]), as shown in Ref. 
0). £(i) = tT{S + t(t)S-\n){n\} is the self-energy (or 
memory kernel) describing the transverse-spin dynamics, 
where £(i) = —iPLQe~ lLClt QLP, P is a projector onto a 
product state of HH and nuclear spins, Q = 1 — P, and 
LO = [H, O] for some operator O acting on the total 
Hilbert space of HH and nuclear spins [TT] . It is conve- 
nient to perform a Laplace transform on Eq. ^ , yielding 
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Figure 1. Decoherence rate I/T2 from Eq. (121 as a function 
of the HH Zeeman energy u) n = g^sB + pi Ahh- For L = 
lOnm and a z — 4nm, we estimate N ~ 7.3 x 10 4 . Inset: I/T2 
for fields up to IT (axes in the same units as in the main 
figure) . 



an algebraic equation of the form 



S + (s + iui n ) = 



(S + )o 

S + ffi(s + lUJ-n 



(7) 



in the frame rotating with frequency u> n . Eqs. |6]) 
and (|7|) are exact equations describing, in general, non- 
Markovian dynamics of the transverse HH-spin compo- 
nent. The structure of the self-energy E(s) is, how- 
ever, very complex, so we have to resort to an approx- 
imation scheme. The energy scales associated with the 
transverse coupling V = (h + S~ + h~S + )/2 are much 
smaller than those associated with the longitudinal cou- 
pling H = (b + h z )S z (see above), and we expand the 
self-energy in powers of hole-nuclear-spin flip-flop pro- 
cesses induced by V: E(s) = E< 2 )(s) + £ (4) (s) + 0(V 6 ). 
Odd orders in V vanish because of the Zeeman mismatch 
between HH and nuclear spins which energetically for- 
bids such processes. For a nuclear spin I of order unity, 
the smallness parameter which controls this expansion is 
given approximately by A± /u) n (see Appendix A of Rcf . 
HB), where A x = y/A 2 CB +A 2 LK . 

We evaluate the second- and fourth-order self-energy 
contributions explicitly, following the procedure de- 
scribed in Ref. [53]. We find, for a homonuclear system 
in the frame rotating with frequency uj n , 



E< a >(» + iw n ) 
E (4) (s-Hw n ) 



4w„ 
. c + c_ 
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(8) 
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where the sums run over all nuclear sites k. We have in- 
troduced c± = 1(1+1) — ((to(to±1))), where I is the nu- 
clear spin, m — —I, and the double angle bracket 
indicates averaging over the 1% eigenvalues m 



We emphasize that the structure of the self-energies 
E^ 2 ) and E^ 4 ) bears some similarity with previous results 
on electron-spin decoherence |26j . However, there are 
two important differences compared to the electron case: 
(i) The appearance of different coupling constants Af. 
and A k in Eqs. ^ and (Job is due to the anisotropy of 
the hyperfine Hamiltonian (I5J) and provides an additional 
smallness factor A±/A z <C 1 (A z = |Ahh|) to the self- 
energy ( 11 ). (ii) The spatial dependence of the transverse 



coupling constants differs from the longitudinal ones due 
to the appearance of excited-state envelope functions. In 
particular, this means that nuclear spins at the edge of 
the QD (rather than in its center as in the electron case) 
couple most strongly to the HH along the transverse di- 
rection - an effect which manifests itself directly in the 
appearance of a distinct minimum in the decoherence rate 
1/T 2 (see Fig. [j). 

We now evaluate the second- and fourth-order self- 
energy in the continuum limit (changing sums to inte- 
grals, see Appendix [5]), following Ref. [25]. Since a z <C L 
(see above), we can perform a two-dimensional limit by 
averaging over the z-dependence in the hyperfine cou- 
pling constants Af, and A^. From Eq. Q, we see that 
the second-order self-energy E^ 2 ) is purely real, leading to 
no decay but a frequency shift Acu = — ReE^ 2 '(s + iu) n ), 



c_ A\ 



167V 



(10) 



where N is the number of nuclear spins enclosed by the 
envelope function. The fourth-order self-energy becomes 



E^ 



.c+c- A±_ A\ 



AN 
1 f i 
dx / dy 
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x (log x) 2 y (logy) 2 
s + i(x — y) 



(11) 



in the continuum limit, where x — exp{— r 2 }, y — 
cxp{— r§}, and = \J x 2 + yf/L (i = 1,2). Here, we 



I^.cbI 2 - 



A fc, LH l 2 since the 



have approximated \A^\ 2 
overlap term vanishes under spatial averaging. The ap- 
pearance of polynomial prefactors r 4 , represented by the 
log functions in the numerator of Eq. (11), is a direct 



consequence of the excited-state envelope functions de- 
scribing the distribution of transverse coupling constants 
Ar within the quantum dot. 

The transverse-spin dynamics of the HH are described 
by the non-analytic structure of the right-hand side of 
Eq. (see, e.g., Ref. [26]). Inserting E< 2 ) and E^ 
into Eq. Q, we find one pole at s ~ iAlu ~ T, whose 
negative real part gives the HH decoherence rate T = 
1/T 2 ~ -ImEW(iw„ + iAw-0 + ) [26J, where 0+ denotes 



a positive infinitesimal. Evaluating Eq. ( 11 ), we find (see 
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Figure 2. Maximum of the decoherence rate I/T2 as a func- 
tion of the QD height a z . For increasing a z , the maximal 
value of I/T2 decreases (see main plot) and the position of 
the maximum is shifted (see inset). 



interaction can be expected to be the dominant source 
of hole-spin decoherence. We emphasize that Eq. ( 12 ) is 
still valid at B = 0, as long as e < 1. For e > 1, Aw ex- 
ceeds the bandwidth of excitations A z /N in the nuclear 
bath, and the perturbation expansion breaks down [26) . 

The degree of band hybridization, and therefore the 
decoherence rate 1/T 2 , depends on the geometry of the 
QD, i.e., on L and a z . For flat QDs the amount of LH 
admixture to the HH states ([I]) is decreased, leading to 
smaller non-Ising terms in the Hamiltonian (p)]). On the 
other hand, the envelope wavefunction of a flat dot en- 
closes less nuclear spins (for fixed L). These two effects 
lead to an increase of the maximal decoherence rate for 
smaller a z (see Fig. |2|, and a shift of its position as a 
function of u) n (see inset of Fig. [2]) . 

We acknowledge discussions with B. Braunecker, W. A. 
Coish, F. Pedrocchi, D. Stepanenko, and M. Trif, and 
funding from the Swiss NSF, NCCR Nanoscience, and 
DARPA Quest. 
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Appendix A: Band hybridization 



1 _ 7TC + C_ Aj_ A\ 

Y 2 AN ~T X ~J£ 

r.1 



x / dxx[\ogx] 2 (x + e)[\og(x + e)} 2 , (12) 



where e = JV|Ao;/Ahh|- The integral in Eq. (12) can 
now be evaluated numerically for any value of e. 

The Zeeman energy of the HH is given by co n = 
ghHsB + pIAna, where —1 < p < 1 is the degree of 
nuclear-spin polarization (along the positive z-direction). 
In Fig. [l] we show the hole-spin decoherence rate I/T2 
as a function of uj n . The non-monotonic behavior of I/T2 
for small uj n appears when e oc l/ui n approaches unity. 

For electrons, a non-monotonic behavior of I/T2 has 
been predicted as well [33] , albeit with a different depen- 
dence on e and around magnetic fields of several Tesla. In 
contrast, for holes, the non-monotonicity occurs at much 
lower fields (B ~ O.lmT for the parameters used in Fig. 
[l] assuming p = 0) , and the rate 1 /T 2 features an addi- 
tional dip which is a footprint of the excited-state enve- 
lope functions appearing in Eq. 0. The huge difference 
in energy scales has very important consequences for the 
tunability of the hole-spin decoherence rate: by increas- 
ing the externally applied magnetic field (or the degree of 
nuclear-spin polarization), it is possible to decrease 1/T 2 
over many orders of magnitude within the experimen- 
tally accessible range of magnetic fields (see inset of Fig. 
[j}. This means that this system offers the possibility 
to entirely 'turn off' hyperfine-associated spin decoher- 
ence. As a consequence, hole-spin quantum dots may be 
operated in a regime where other interactions, such as 
spin-orbit or direct nuclear dipole interactions, will be 
the dominant source of spin decoherence and will there- 
fore become experimentally observable. On the other 
hand, for small uj n , the hybridization- induced transverse 



Here we show how to derive the hybridized heavy-hole 
states given in Eq. (1) of the main text. We start from 
the bulk version of the Kane Hamiltonian, which reads 
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where the relevant blocks are given by 
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B = -e[( 7 i ^ 
E = Pk+, 
F = e[ l2 {k 2 x - 
J = ej 3 k-k z 



72)(^ + ^) + ( 7l -2 72 )fc 



2 1 

z J - 



2i^k x k 



(A3) 



Here, e = fi 2 /2rao, mo is the free-electron mass, 
k± = k x ± iky, 7j denote the Luttinger parameters, and 
P is the inter-band momentum. -Hcb, -Hlh, and -ffso de- 
scribe unperturbed electron states in the CB, LH and SO 
bands, while V 2 , V3, and V 6 describe the CB-LH, CB-SO, 
and LH-SO coupling, respectively, and thus do not con- 
tribute to the first order of the perturbation expansion 
carried out in Eq. (A8) below. We have neglected terms 
that are proportional to C, B 7v , and B^ v in Winkler's 
notation [23] since, due to their smallness, they will not 
lead to considerable corrections for our purposes. 
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We choose a parabolic confinement potential 



Vdot = V z (z) + V xy (x,y), 



V z (z) 



mj_uj | 2 

zr^Z , 



K?/(a;,y) 



(x 2 +y 2 ), (A4) 



where u;j_ = h/(m±a 2 ) and oj| = h/(m\\L 2 ) with HH 
effective masses toj^ = to /(7i — 2j 2 ) and my = Too/ (71 + 
72) along the growth direction and in the plane of the 
dot, respectively, and a z and L denote the corresponding 
confinement lengths. 

In order to calculate the hybridized HH states per- 
turbatively, we take the ground-state eigenket |0h° h ) of 
ffnH + Vdot as the unperturbed envelope function, whose 
position representation is given by 



with harmonic-oscillator eigenfunctions 



1 
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2 ~A 
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(A6) 
(A7) 



Recall that the position representation of the envelope 
functions in band a is given by (v\<p 1 ^) — 4> l ^-{z)4& (x, y) 
(i = 0,1, 3 = 0,±), where <t$(x,y) = 0° 
fi ll (x,v) = (<Pi(x)f a (y)±i<Pl(x)^(y))/V2, ajid 4>»(x) 
is the n th harmonic-oscillator cigenfunction in band a. 
The hybridized HH states can now be evaluated by a 
perturbation expansion (in inverse energy splitting of 
the bands) as follows: defining a two-spinor I^hh) = 
(I^hh)' I^hh))*' the hybridized HH spin states are the 
components of 
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so 



l^soX^ol^ Whh), 



where I, to, and n label the electronic states in the CB, 
LH, and SO band, respectively, E g is the fundamental 
band gap, A LH is the HH-LH splitting (see Eq. ( |A14[ )), 
and Ago is the spin-orbit-induccd LH-SO band splitting. 

Due to the cylindrical symmetry of the dot, only one 
matrix element per band yields a non-vanishing contri- 
bution to the sum. For instance, (^c°bI^±I^h°h) = due 
to the vanishing angular integral, but (0 C g|fc± |</>hh) 7^ 0' 



where (r| 



O^OTfov) with 



x ±iy 



4>CB{x,y). 
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L = La I + 72 is the renormalized lateral confine- 
ment length of a CB electron (due to the difference in ef- 
fective masses between particles in the CB and HH band) , 



where a = \J mo / m* with the free-electron mass Too and 
the effective CB electron mass m* . The excited-state 



envelope function ( A9 ) features a polynomial prefactor 



which gives rise to the appearance of logarithmic terms 
in the self-energy ( |B3[ ) and, eventually, to a minimum of 
the decoherence rate 1 /T2 as a function of uj n (see Fig. 1 
of the main text). The matrix element 
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-72 J 



272) 

(A10) 

differs from unity because of the different effective CB 
and HH masses which lead to different spring constants 
of the harmonic-oscillator potential in the CB and HH 
bands. Similarly, for the LH contribution, the only non- 
vanishing term in the sum in Eq. ( A8 1 is given by 



^lhI^±^zI^h°h)> an d the overlap integral reads 



Ah = <0 H hIO = 
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The HH-LH splitting Alh is given by the difference in 
ground-state energies in the HH and LH bands: Alh = 
£hh - s lh, where 
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This leads to a HH-LH splitting of 
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Carrying out the perturbation expansion following Eq. 



(A8), inserting the quantities given above, and neglecting 



the small SO band contributions induced by V5, we arrive 
at the hybridized states given in Eq. (1) of the main 
text. Note that the Kane Hamiltonian ( Al ) is written in 



the basis of CB, HH, LH, and SO Bloch states, so that 
these Bloch states only appear implicitly in Eq. (A8|, 
but explicitly in the final form given in Eq. (1) in the 
main text. In order to produce explicit numbers, the 
Bloch amplitudes have been approximated by a linear 
combination of n — 4 hydrogen-like atomic orbitals, such 
as in Ref. [5]. 

We have also neglected HH-LH coupling terms induced 
by the F matrix element in V4 since such terms are para- 
metrically suppressed by a factor of a z /L <C 1 with re- 
spect to the couplings induced by the J matrix element. 
Note that the LH (pseudospin) Bloch functions near the 
r-point, 



|«LH±)|± 



LH/ 



(V2\ Pz )\t,l)T\p±)\U))/V3, (A15) 



already include (real) spin-up and spin-down states, such 
that the J matrix element, although diagonal in pseu- 
dospin space, induces off-diagonal couplings in the effec- 
tive Hamiltonian (Eq. (5) of the main text). Therefore, 



6 



0.4 



> 

^ 0.2 



s 



. ^CB 

- - A hH 


/ 

s 

* 

* 

* 

* 

* 










i i 


1 1 1 



a z [nm] 

Figure 3. Conduction-band (CB) and light-hole (LH) contri- 
butions to the non-Ising hyperfine coupling A± as a function 
of the perpendicular dot size a z . For this plot, we have chosen 
a GaAs dot of lateral size L = lOnm. 



We now replace sums by integrals according to vq Y2k ~~ * 
/ d 3 r, where vq is the volume occupied by one nucleus 
and the integration is carried out over all space. Aver- 
aging over the z-dependence of the coupling constants, 
A(x,y) = dzA(x,y,z) we write X^ 4 - 1 in terms of a 
four-dimensional spatial integral: 



£ (4) (s + iw„) ~ -i 



c+c_ A 4 



271 



d6 l 



4ir 2 N lu 2 A : 
d9 2 



dr 1 r 1 / dr 2 r 2 

n JO 



o 



s + i{e r i — e 2r 2 ) ' 



(B2) 



where we have introduced spherical coordinates Xi = 
ri cos 9i, yi = r i sin 8i, and where s = sN/A z . We have 
also rewritten Nvq = ira z L 2 (the volume of the QD). 
The angular integrals simply contribute a prefactor of 
Ait 2 . The radial integrals can be solved by introducing 
new variables x = e~ Tl , y = e~ V2 , such that 



the F matrix element does not introduce qualitatively 
new contributions to the effective Hamiltonian. 

The degree of band hybridization is described by the 
dimensionless prefactors 



Acb = 



ijhBP 
\/2LE n 



2V2 72 (I 2 -ai; 



(A16) 



We see that these prefactors, and therefore also the hy- 
perfine coupling strengths Aqb oc |Acb| 2 and Alh oc 
|Alh| 2 j have a strong dependence on the QD geometry 
(i.e., on L and a z , therefore the non-Ising terms in Eq. 
(5) of the main text can, in principle, be 'tuned' to be 
more CB or LH-like (see Fig. [3]). It should be empha- 
sized, however, that CB and LH admixture lead to qual- 
itatively similar contributions to the effective Hamilto- 
nian, Eq. (5) in the main text. For a quantum dot with 
L — lOnm and a z — 2nm, we estimate |Acb| — 0.04, 
|A LH | - 0.11 for GaAs. 



Appendix B: Continuum limit 

Here we show in detail how to calculate the decoher- 
ence rate Y — I/T2 from the fourth-order self-energy in 
the continuum limit. We start from Eq. (9) of the main 
text, 



S (4) (s + iw„) ~ -i 



C_|_C_ 

4w 2 



Ki 2 i4 12 



E 

ki,k 2 



Ai 



(Bl) 



£ (4) (s-Hw„) ~ -i 



c+c^ A\ 



x / dx 

'0 ./ii 



1 x (log x) 2 y (logy) 2 
dy 



s + i(x — y) 



(B3) 



The right-hand side of the equation of motion in 
Laplace space (Eq. (7) of the main text), 



S + (s + iui n ) = 



(S + )o 



s + iS(s + iuj n ) 



(B4) 



features a pole at s — iAuj — Y [HI [55], where Y = 
— Im£( 4 )(/u/„ + iAcj — + ) and where + denotes a posi- 
tive infinitesimal. Evaluating E^ 4 ' at s = iu n + iAw — + 
and using 



(B5) 



where V denotes that the principle value should be taken 
in any integration involving the above expression, we ar- 
rive at an integral of the form 



dx \ dyf(x,y) 
n Jq 

x \V 



1 



m5(x -y + Aw) . (B6) 



x — y + Aw 



Taking the imaginary part according to Y = 
-Im£W(iu/„ + iAuj - 0+) leads directly to Eq. (12) of 
the main text. 
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